Introduction to Data

Learning objectives

By the end of the lab, you will be able to …

  • implement basic variable manipulation
  • create useful frequency tables

Code-along 02

Download and open code-along-02.qmd

Intro. to data

Packages

Load the standard packages.

library(here)
library(tidyverse)
library(gssr)
library(gssrdoc)


Install and load the summarytools package.

install.packages("summarytools")


library(summarytools)

Load your data & codebook

# Get the data only for the 2024 survey respondents
gss24 <- gss_get_yr(2024)

# Load the codebook
data(gss_dict)

Operators in R

Operators in R are symbols directing R to perform various kinds mathematical, logical, and decision operations. A few of the key ones to know before we get started:

To test equality or inequality:
==, !=, >, >=, <, <=

To indicate “and”, “or”, and “not”:
& | !

Assigning values to various data objects: <- -> =

function(argument)

Functions are (most often) verbs, followed by what they will be applied to in parentheses:


do_this(to_this)
do_that(to_this, to_that, with_those)

Variables

Remember, you can access the variables (i.e., columns) using the $ operator, as shown using the table() function.


The variable names are case sensitive. In this dataset, all variables are lowercase.

table(gss24$fefam)

  1   2   3   4 
167 492 899 604

195 respondents were coded as 2 on this variable. What does that mean?

Variable types

Frequency Distributions

Variable Descriptions

Let’s familiarize ourselves with the premarsx and polviews variables.


In the console, type ?premarsx and hit enter. The Help pane will show you the question text, response options and values.


Now, do the same for polviews.

Table

Run this code to see the frequency table for the premarsx variable. Then, add a line below to also see a table for the polviews variable.

table(gss24$premarsx)

   1    2    3    4 
 357  122  258 1378 


table(gss24$polviews)

   1    2    3    4    5    6    7 
 140  421  368 1148  381  516  186

Cross-tabs

The table command also let’s you create a table with two variables.

# 1st variable is the rows, 2nd variable is the columns.
table(gss24$premarsx, gss24$polviews)
   
      1   2   3   4   5   6   7
  1   8  11  13  78  45 132  52
  2   1  10  10  44  18  25   8
  3   3  26  29  91  41  43  16
  4  91 227 187 488 148 145  38

Labels

Use haven::as_factor to see the value labels instead of the value numbers. Then, do the same for polviews.

table(as_factor(gss24$premarsx))

                 always wrong           almost always wrong 
                          357                           122 
         wrong only sometimes              not wrong at all 
                          258                          1378 
                        other                           iap 
                            0                          1126 
                   don't know            I don't have a job 
                           50                             0 
                  dk, na, iap                     no answer 
                            0                             6 
                not imputable                       refused 
                            0                             0 
               skipped on web                    uncodeable 
                           12                             0 
not available in this release    not available in this year 
                            0                             0 
                 see codebook 
                            0

Labels

table(as_factor(gss24$polviews))

            extremely liberal                       liberal 
                          140                           421 
             slightly liberal  moderate, middle of the road 
                          368                          1148 
        slightly conservative                  conservative 
                          381                           516 
       extremely conservative                    don't know 
                          186                            99 
                          iap            I don't have a job 
                            0                             0 
                  dk, na, iap                     no answer 
                            0                            20 
                not imputable                       refused 
                            0                             0 
               skipped on web                    uncodeable 
                           30                             0 
not available in this release    not available in this year 
                            0                             0 
                 see codebook 
                            0

Better Labels

Let’s clean up the levels for premarsx.

gss24$premarsx <- zap_missing(gss24$premarsx)
gss24$premarsx <- as_factor(gss24$premarsx)
table(gss24$premarsx)
1
Get rid of all the ‘missing’ levels (just missing)
2
Apply the labels instead of numeric values 

        always wrong  almost always wrong wrong only sometimes 
                 357                  122                  258 
    not wrong at all                other 
                1378                    0

Better Labels cont.

Let’s get rid of the empty levels in premarsx.

gss24$premarsx <- droplevels(gss24$premarsx)
table(gss24$premarsx)

        always wrong  almost always wrong wrong only sometimes 
                 357                  122                  258 
    not wrong at all 
                1378

Manipulating Variables

For polviews, let’s combine categories to ease interpretation. This is easiest when the levels are numeric.

Let’s remind ourselves what the values correspond with each label.

table(as_factor(gss24$polviews, levels = "both")) # both shows value and label

             [1] extremely liberal                        [2] liberal 
                               140                                421 
              [3] slightly liberal   [4] moderate, middle of the road 
                               368                               1148 
         [5] slightly conservative                   [6] conservative 
                               381                                516 
        [7] extremely conservative                    [NA] don't know 
                               186                                 99 
                          [NA] iap            [NA] I don't have a job 
                                 0                                  0 
                  [NA] dk, na, iap                     [NA] no answer 
                                 0                                 20 
                [NA] not imputable                       [NA] refused 
                                 0                                  0 
               [NA] skipped on web                    [NA] uncodeable 
                                30                                  0 
[NA] not available in this release    [NA] not available in this year 
                                 0                                  0 
                 [NA] see codebook 
                                 0

Manipulating Variables

gss24 <- gss24 |>
 mutate(pol3cat = case_when(
   polviews >= 1 & polviews <= 3 ~ "Liberal",
   polviews == 4 ~ "Moderate",
   polviews >= 5 & polviews <= 7 ~ "Conservative",
   TRUE ~ NA_character_),
  pol3cat = factor(pol3cat,
                 levels = c("Liberal", "Moderate", "Conservative"))
  )
1
Save over the dataset with an added variable.
2
Creates a new variable by assigning labels based on values of polviews
3
Set everything else to “missing”
4
Convert the variable to factor with specified level order


Pipes

can be written as |> or %>%

Frequency Table

Always double check your work.


table(gss24$polviews, gss24$pol3cat)
   
    Liberal Moderate Conservative
  1     140        0            0
  2     421        0            0
  3     368        0            0
  4       0     1148            0
  5       0        0          381
  6       0        0          516
  7       0        0          186

Relative Frequency Table

Make a frequency table. One of summarytools main purposes is to help cleaning and preparing data for further analysis. Pay attention to the missing values. Then, do the same for premarsx.


freq(gss24$pol3cat) 
Frequencies  
gss24$pol3cat  
Type: Factor  

                     Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
------------------ ------ --------- -------------- --------- --------------
           Liberal    929     29.40          29.40     28.07          28.07
          Moderate   1148     36.33          65.73     34.69          62.77
      Conservative   1083     34.27         100.00     32.73          95.50
              <NA>    149                               4.50         100.00
             Total   3309    100.00         100.00    100.00         100.00

Relative Frequency Table

freq(gss24$premarsx) 
Frequencies  
gss24$premarsx  
Type: Factor  

                             Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
-------------------------- ------ --------- -------------- --------- --------------
              always wrong    357     16.88          16.88     10.79          10.79
       almost always wrong    122      5.77          22.65      3.69          14.48
      wrong only sometimes    258     12.20          34.85      7.80          22.27
          not wrong at all   1378     65.15         100.00     41.64          63.92
                      <NA>   1194                              36.08         100.00
                     Total   3309    100.00         100.00    100.00         100.00

Pretty Tables

Using report.nas = FALSE suppresses the missing data.
The headings = FALSE parameter suppresses the heading section. Do the same for premarsx.


freq(gss24$pol3cat, report.nas = FALSE, headings = FALSE) 

                     Freq        %   % Cum.
------------------ ------ -------- --------
           Liberal    929    29.40    29.40
          Moderate   1148    36.33    65.73
      Conservative   1083    34.27   100.00
             Total   3160   100.00   100.00

Pretty Tables

freq(gss24$premarsx, report.nas = FALSE, headings = FALSE) 

                             Freq        %   % Cum.
-------------------------- ------ -------- --------
              always wrong    357    16.88    16.88
       almost always wrong    122     5.77    22.65
      wrong only sometimes    258    12.20    34.85
          not wrong at all   1378    65.15   100.00
                     Total   2115   100.00   100.00

Cross-tab

The table() function gives us the frequencies.


table(gss24$premarsx, gss24$pol3cat)
                      
                       Liberal Moderate Conservative
  always wrong              32       78          229
  almost always wrong       21       44           51
  wrong only sometimes      58       91          100
  not wrong at all         505      488          331


We want to add the column percentages…

Relative Frequency Cross-tab

ctable(gss24$premarsx, gss24$pol3cat,
       prop = "c",
       format = "p",
       useNA = "no")
1
Change from table() to ctable().
2
The “c” gives column %; “r” would give row %.
3
This adds the % symbols to the table.
4
Exclude the missing levels from the table.